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Abstract 

We study the baby Skyrme model in a pure AdS background without a mass 
term. The tail decays and scalings of massless radial solutions are demonstrated to 
take a similar form to those of the massive flat space model, with the AdS curvature 
playing a similar role to the flat space pion mass. We also numerically find minimal 
energy solutions for a range of higher topological charges and find that they form 
concentric ring-like solutions. Popcorn transitions (named in analogy with studies 
of toy models of holographic QCD) from an n layer to an n+l-layer configuration are 
observed at topological charges 9 and 27 and further popcorn transitions for higher 
charges are predicted. Finally, a point-particle approximation for the model is 
derived and used to successfully predict the ring structures and popcorn transitions 
for higher charge solitons. 
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1 Introduction 


The Skyrme model [lj is a non-linear theory of pions in (3 + 1) dimensions, admitting 
soliton solutions called Skyrmions. It has been derived as a low-energy effective field 
theory of QCD in the large colour limit, and has also been used in holographic models 
such as the Sakai-Sugimoto model where Yang-Mills Chern-Simons instantons in a 
(4+ l)-dimensional bulk spacetime are dual to (extended) Skyrmions on the boundary. A 
key feature of these models is that the bulk spacetime is AdS-like, possessing a conformal 
boundary and a finite, negative scalar curvature. 

Solitons in pure Anti de-Sitter spacetimes are also of interest. It has been shown that 
Skyrmions with massless pions in hyperbolic space are closely related to Skyrmions with 
massive pions in Euclidean space HE]- Since constant time slices of AdS spacetimes 
are hyperbolic we may expect similar results in pure AdS. In addition, monopoles and 
monopole walls have been studied in AdS HE], motivated as a magnetic version of 
holographic superconductors. 

The baby Skyrme model |£] is a (2 + 1) dimensional analogue of the Skyrme model. 
Baby Skyrmions have recently been used to study low-dimensional models of the Sakai- 
Sugimoto model in the context of dense QCD [9] HU] . In these toy models a series of 
phase-transitions were observed in which infinite chains of solitons split into multiple 
layers with increasing density. These were dubbed popcorn transitions, and the extra 
layers were found to be separated in the holographic direction. 

Here we investigate the baby Skyrme model in a pure AdS background and study the 
resulting soliton and multi-soliton solutions. The low dimensionality of the model makes 
full numerical field computations viable, and we find that the curvature of the spacetime 
allows us to find soliton solutions even without a pion mass term. Multi-solitons beyond 
topological charge 5 = 3 are found to take the form of ring-like structures, with popcorn- 
like phase transitions to multi-layered rings occurring as the topological charge increases. 
Inspired by methods used to study aloof baby Skyrmions HU we will derive a point- 
particle approximation to study these phase transitions for higher topological charges. 


2 AdS spacetime in (2 + 1) dimensions 

AdS is the maximally symmetric spacetime with Minkowskian signature and constant 
negative curvature, and in this paper we will be interested in the (2 + l)-dimensional 
case. The metric, given in sausage coordinates, can be written 

ds2 = - ( i ^ y ) dt% + (1 ^2)2 ( dx2 + d V 2 ) • (!) 

Here, L is the AdS radius (related to the cosmological constant A via A = —1/L 2 ) and 
r = \Jx 1 + y 2 G [0,1) is a radial coordinate. These coordinates are useful numerically 
as they are global coordinates over a finite range. Later we will use that the geodesic 
distance between two points x and y in this spacetime is given by 

d (x,y) = L coslr' (l + . (2) 


2 






The following global coordinates give another useful way of writing the AdS metric: 

ds 2 = — cosh 2 y- dt 2 + dp 2 + L 2 sinh 2 y dO 2 , (3) 

±j ±J 


where p > 0. These coordinates are useful since the radial coordinate p coincides with 
the geodesic distance from the origin in this model. The global and sausage coordinates 
are related by r = tanh y-. 

Since constant time slices of AdS are hyperbolic space, it will be useful to define a 
hyperbolic translation in sausage coordinates. A translation sending the origin to some 
point a is given by 


(1- 

a 

\ 2 )x + (1 + 2x 

• CL + 

\x\ 

2 )a 

1 -\- 2 jX • CL -\- 

a 

2 

\x\ 

2 


(4) 


It should be noted that this is an isometry of only the constant time slices of the spacetime, 
rather than the full spacetime, since the time component of the metric is dependent on 
r. 

As a final note, the Ricci scalar curvature of AdS can be calculated as R = —6 /L 2 . 
In the limit L —>■ oo this curvature vanishes and we recover flat space. 


3 The AdS baby Skyrme model 


We will be interested in studying soliton solutions of the baby Skyrme model, given by 
the action 

s B3 = -\f (d,4> • 04 + ^( 34 x 84 ) • (84 X 84) 

+ 2m 2 (l — <f> • n)^ a/— det ~g dx dy dt . (5) 


The field 0 = (0i, 02, 03) is a three component unit vector. The first term is that of the 
0(3)-sigma model, and the second term is the baby Skyrme term with constant coefficient 
k 2 . The third term is a potential term containing the pion mass parameter rn (named in 
analogy with the pion mass term in the full Skyrme model), and a constant unit vector 
n. Greek indices run over spacetime coordinates t, x and y, and later we will use Latin 
indices to run over the purely spatial coordinates. The symmetry of the action is 0(3) 
for m = 0 and 0(2) for m/0. 

The associated static energy of this model is 


E B S = \f 1 T ^{\dA\ 2 + \d y 4?* 


which yields the equations of motion 


2^,2 


k 2 (1 — r 2 ) 
4L2“ 


■\d x <f> x d y 4 >| 2 

8L 2 m 2 w i i 

+ (1 _ r 2 ) 2 ( 1 ~^- n ))^^^ 


( 6 ) 


1 / i~2 ( i _ y,2\2 

diji = (% <j --— (0 x di<t> +- — dj<j){<t> ■ 84 x 84) 


1 — r 2 
= m 2 \f— det ~gn x 0. 


(7) 
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For finite energy we then require <j> —>■ n as r —» 1. Without loss of generality we 
can choose n = (0, 0,1), and in the massless (m = 0) case this choice of boundary value 
breaks the 0(3) symmetry of the model to 0(2). We can identify points on the boundary 
r = 1 and treat the pion fields as maps (p : S 2 —» S 2 , giving rise to an associated winding 
number and topological charge 

-B = ; J<t> ■ (d x <f> x d y (t>) dx dy (8) 

which we identify with the baryon number of the configuration. By noting the inequalities 

\d x (j)±<l> x dy(f>\ 2 > 0, ( 9 ) 

we obtain the Bogomolny bound Ebs > 4tt | i? |. 

4 Radial baby Skyrmions in AdS 

We begin by discussing some properties of radially symmetric solitons in our model, 
working in global coordinates 0- Due to the principle of symmetric criticality and the 
symmetries of both AdS and the action, we would expect static B = 1 solitons in our 
model to posses radial symmetry and be centred at the origin p = 0. Static, radially 
symmetric configurations with topological charge B are given by the hedgehog ansatz 

(f> = (sin/(p) cos {B{9 - ip)), sin/(p) sin (B(9 - ip)), cos f (p)) (10) 

where f(p) is some profile function satisfying /(0) = n, /(oo) = 0. ip is some constant 
internal phase which we can set to zero here due to symmetry, although internal phase 
differences will become important when we consider multi-solitons. Substituting into 
and performing the coordinate transformation yields the static energy 

Erad = ~Yj Q sinh y ^/ /2 + ^ + ^ 2 / /2 ) + 2 ™ 2 (1 ~ cos /) j d P • ( n ) 

We can numerically find the profile functions f(p) for different values of B by min¬ 
imising © using a modified gradient flow method. 

In other baby Skyrme models it has been found that radial solitons can be well- 
approximated by flat-space instantons of the 0(3)-sigma model. This approximation 
has enabled an investigation of how the soliton sizes p scale with the baby Skyrme 
parameter k. Unfortunately 0(3)-sigma instantons in AdS do not have a convenient 
closed form expression, due to the presence of the non-constant time component of the 
metric, preventing us from analytically exploring this relationship. We can nevertheless 
perform a numerical investigation; since the profile function interpolates between 7r and 
0 we can define the size of the soliton as p : f(p) = 7 t/ 2 . 

Using this definition it is straightforward to numerically find how p scales with k and 
L for different values of B. In the massless case we find the leading-order dependence is 
p ~ VkL for small k/L; looking at © we see that nonlinear effects will dominate when 
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k/L is large. Comparing this to the scaling of baby Skyrmions with mass parameter 
m > 0 in flat space, p ~ we can see that the curvature of AdS space can be 

interpreted as adding an effective pion mass. 

Finally, we can calculate the leading-order decay of the soliton tails near the boundary 
of our space. Assuming a radial ansatz, and using the fact that f(p ) —>■ 0 as p —>■ oo, we 
can linearise the equations of motion to obtain 


L 2 sinh — /" + 2 L cosh —/' 
L J L J 




+ m L 


2 t 2 


f = o. 


( 12 ) 


In the limit p —>■ oo we can obtain the asymptotic tail decay as 

f(p) ~ e -(!+FT WZ*) P /L ^ ( 13 ) 


independently of k or B. This relation also holds in the special case m — 0. It is 
interesting to note that, unlike baby Skyrme models in flat space, the addition of a mass 
term is not required for the soliton tail to decay exponentially. This is in contrast with 
baby Skyrmions in flat space which have large-radius asymptotic tail decays 


f(r) ~ 



if m — 0 
if m 7 ^ 0 . 


(14) 


In fact, these flat-space tail decays can be obtained from the linearised asymptotic equa¬ 
tions of motion by carefully taking the limit L —y oo, as expected. 


5 Multi-solitons in AdS 

As we shall see, for topological charges B > 3, radial baby Skyrmions in AdS no longer 
provide global minima for the energy functional ([ 6 ]). Analytic investigation of higher 
charge solitons is difficult, but we are able to perform full numerical minimisations to 
seek out local and global energy minima. The numerical results in this section were ob¬ 
tained by performing a modified gradient flow method on the energy (| 6 ]) with parameter 
values k = 0.1, L = 1, m = 0 on a grid with 501 x 501 gridpoints in sausage coordinates. 
Derivatives were calculated used fourth-order finite difference approximations. In addi¬ 
tion, the numerical results were verified by applying a fourth-order Runge-Kutta method 
to the full dynamical equations of motion. 

In order to find local energy minima we investigated a range of different initial condi¬ 
tions for our minimisation algorithms. Initial conditions were primarily generated using 
the product ansatz: if we project a field (f> onto the Riemann sphere by 

W = + i<h 

l-<h ’ 

then we can generate a field configuration composed from two fields Wi, IF? by writing 

WiW 2 

= 

Wi + W 2 



(15) 
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Figure 1: Plots of the numerically calculated AdS baby Skyrmions for topological charges 
l<£<20forK = 0.1,L = l and m = 0. We plot <p 3 in sausage coordinates as a visual 
representation of the solitons. 
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Positions of solitons can be identified with points where <f> = (0,0,—1) or equivalently 
W = 0, and W = oo whenever (f> = (0,0,1), the boundary value. Our composed field 
has zeroes at the zeroes of W i and W ?, so we see that the product ansatz gives us a 
field with the correct topological properties as the superposition of two individual fields. 
Combining this with the hyperbolic translations Q allows us to generate initial conditions 
that resemble collections of radial baby Skyrmions placed in different positions on our 
grid. We also used perturbed radial fields as initial conditions. 

Figure [l] shows colour contour plots of 0?, for numerically found global energy minima 
for topological charges 1 < B < 20, with energies given in Table [lj For the pion mass 
term used in this paper, radial solutions are preferred up to charge B = 3, in contrast to 
baby Skyrmions in flat space which only have radial energy minima for B < 2. 

For 4 < B < 7 the baby Skyrmions form regular polygons, with soliton positions at 
the vertices. The relative phase differences between neighbouring soliton positions is 7r 
or 7r ± it/B for even or odd charges respectively. 

At charge B = 8 the ring structure deforms due to the centralising force induced 
by the AdS metric. For B > 9 the solutions form multi-layered concentric rings, and 
the central layers for 9 < B < 16 resemble (normally slightly deformed) radial solutions. 
These deformations appear to match the symmetries of the outer rings. We denote multi¬ 
layered ring structures as {ni, ri 2 , ,... }, where n t denotes the topological charge in the 
ith ring, counting out from the origin. 

The transition from a single ring to a multi-ring structure with increasing baryon 
number is reminiscient of the popcorn transition observed in low-dimensional analogues 
of the Sakai-Sugimoto [9, l()j model of holographic QCD, in which finite density chains of 
solitons pop out into a holographic direction at some critical baryon density. In light of 
this, it would be interesting to see if the AdS baby Skyrme model studied here possesses 
further popcorn-like transitions beyond B = 20. 

This is a potentially difficult task to accomplish numerically. As higher charge so¬ 
lutions are investigated the number of local energy minima increases drastically. This 
requires a very large number of initial conditions to be tested in order to gain confidence 
that a global energy minimum has indeed been found. 

In the following sections we will formulate an approximation to our system in which 


charge B 

E/AttB 

form 

charge B 

E/AnB 

form 

i 

1.2548 

{i} 

n 

1.5368 

{2, 9} 

2 

1.2312 

{2} 

12 

1.5554 

{2, 10} 

3 

1.2878 

{3} 

13 

1.5788 

{2, 11} 

4 

1.3384 

{4} 

14 

1.6017 

{2, 12} 

5 

1.3725 

{5} 

15 

1.6250 

{3, 12} 

6 

1.3886 

{6} 

16 

1.6481 

{3, 13} 

7 

1.4263 

{7} 

17 

1.6714 

{4, 13} 

8 

1.4541 

{8} 

18 

1.6914 

{4, 14} 

9 

1.4888 

{1,8} 

19 

1.7107 

{5, 14} 

10 

1.5157 

{1,9} 

20 

1.7276 

{6, 14} 


Table 1: Energies per charge and forms for AdS baby Skyrmions with topological charge 
1 < B < 20. 
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the solitons are modelled by point particles in a gravitational potential with an inter- 
soliton interaction. We will use this model to predict the baryon numbers at which 
further popcorn transitions occur, and use the predictions as a guide to finding further 
global minima in the full model. 


6 AdS baby Skyrmions as point particles 

The numerical results from the previous section are reminiscent of the results of circle 
packings within a circle [12J: finding the minimal area circle within which you can pack 
B congruent circles. Solutions tend to be arranged as rings of points separated by at 
least the diameter of the circles (~ 2/j, for our baby Skyrmions). However, using this 
as an approximation to the AdS baby Skyrme model presents some problems. Firstly, 
the popcorn transitions occur too early, the first two at B = 7 and B = 19. This is 
likely due to the malleable and overlapping nature of the baby Skyrmions. This problem 
still occurs even when the circle packings are formulated in hyperbolic space. However 
it does suggest that a point-particle approximation could be able to qualitatively predict 
the form of solutions, if a better representation of the soliton interactions can be found. 

In order to derive a better point-particle approximation it will be necessary to ob¬ 
tain numerical approximations to the effective gravitational potential and inter-soliton 
interaction. We will make use of the hyperbolic translations @ , and assume that the 
energies of radial fields translated in this way can approximate the energies of constituent 
parts of multi-soliton configurations located at different points on our grid. It should be 
clarified that such fields are not solutions to the equations of motion since the hyperbolic 
translations are not isometries of AdS. 


6.1 The gravitational potential 

We begin by deriving an approximate gravitational potential from the metric ([Tj) . The 
geodesic equations associated with this metric are given by 


t" = - 


x" = — 


1 - r 4 


(xx' + yy')t ', 


X(l+r 2 ) 2 , 2x , a 2 4x 


V 


L 2 ( 1 - r 2 ) 
y(l + r 2 ) , 


0 t'Y + 


1 — r‘ 


-Xv'f - 


1 — r 


-, x v 


(17) 


(‘T+AAm 2 - t%*v. 


L 2 (l _ r 2 ^ 1 — r 2 1 — r 2 

where primes denote differentiation with respect to proper time. In the non-relativistic 
limit y' <C t' we can write 


x 


x't" x h 


x — 


x(l + r 2 ) _ x 


— --- 




(18) 


{t') 2 ( t'Y {t') 2 L 2 ( 1 — r 2 ) r 

where we have implicitly defined the gravitational potential <I>(r). Integrating along with 
the condition $(0) = 0 gives 


d>(r) = f 
J 0 


R{1 + R2) dR = -dx(r 2 + 2 log (1 - r 2 )). 


L 2 (l - R 2 ) 


2 L 2 


(19) 
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Figure 2: Numerical and analytical approximations for the gravitational potential induced 
by the AdS metric. The analytic approximation is <f>(r) = —a(r 2 + 21og(l — r 2 ))/2L 2 
with a = 64.3 and L = 1. Energies are given in units of 47T. 

To fit this potential to the AdS baby Skyrme model we are required to multiply it 
by some constant factor a. We obtain a numerical approximation for the gravitational 
potential by evaluating the energies of translated B = 1 radial fields and subtracting 
off the energy of the true B = 1 solution. We then fit a to this data by performing 
a least-squares fit. We fit only within the radius r = 0.6 because full numerical local 
minima, even for high charges, do not lie much beyond this radius, and the hyperbolic 
translated B = 1 fields become less accurate as approximations near the edge of the disc. 
In addition, we investigated different radii to fit our data to, and found that choosing 
values larger than 0.6 resulted in point-particle approximations that did not successfully 
estimate the full numerical results. 

Figure [2] shows the analytic potential for a = 64.3 compared to the numerical approx¬ 
imations with k = 0.1, L — 1. The curves are in close agreement for radii in the range 
r e [0,0.6], but diverge as r increases further, as expected. 

6.2 The inter-soliton interaction 

We can obtain a numerical approximation for the inter-soliton interaction in a similar 
way: we can numerically calculate the energy of a product ansatz of two translated 
B = 1 solitons, subtract off the potential energies associated with each component soliton 
(according to the numerical approximation above) and the energy of two B = 1 solutions, 
and plot the resulting energy as a function of the geodesic separation of the positions of 
the solitons. 

These static approximations are shown as red curves in Figure [4| where the upper 
curve represents the inter-soliton energy of a pair of solitons in phase, and the lower curve 
represents solitons out of phase. Relative phase differences are calculated with respect 


6 - 


c 

<D 

+-> 

o 

Q. 

ro 

c 
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to the geodesics that the particles lie on (see Fig [3]), so that a pair of soliton field with 
internal phases ip a and 'Wi, have a relative phase difference x — X.(V'a, V'b)- In keeping with 
results found in other spacetimes, we find that pairs of solitons at large separations are in 
the maximally repulsive channel when they are in phase (x = 0), and in the maximallFor 
higher topological chargesy attractive channel when they are out of phase (x = vr). 

With this numerical data as a guide we can fit the out-of-phase interaction using a 
Morse potential of the form 

U„(p) = D (e^-p/pe) _ 2e a(i-p/p e )) } (20) 

where p is the geodesic separation between the solitons, D is the depth of the potential at 
its minimum, p e is the separation at which the potential is minimised, and a is a parameter 
controlling the width of the potential. These parameters can be fit by performing a 
least-squares fit with the numerical data. Since the product ansatz is only valid for well- 
separated solitons, we fit the data in the region where the separation between the solitons 
is greater than twice the size of a single soliton. This yields parameter values D = 0.76, 
p e = 0.73 and a = 1.13. 

In order to introduce a dependence of the potential on the relative phase difference of 
the two solitons, Xi we assume that the solitons are at their most attractive when they 
are out of phase, and their most repulsive when they are in phase i.e. 

U x (p) = D (e^-p/pe) + 2 cos (x)e a(1_p/pe) ) . (21) 

Plots of U( } and U n are given by the upper and lower blue curves in Figure [4] respectively. 

7 Baby Skyrmion rings and shells 

We now seek energy minima of the point-particle approximation above for various values 
of B. For a configuration of B solitons with disc coordinates x a and internal phases xp a 
(where 1 < a < B) we seek to minimise the energy 

Eb = ^ ] |^(bi) + ^ y ^ E x (^ a ^ b ){d{x a ,Xb)) \ , (22) 

a =1 \ b>a / 

where r a = |x a |, <F(r) and U x (p) are the potentials given above and d(x a ,Xb ) is the 
geodesic distance ([2]) between points x a and x^. 

We minimise the point-particle energy using a multi-start stochastic hill-climbing 
method, where a randomly generated initial condition is allowed to relax iteratively. For 
each topological charge we used randomly generated initial conditions, and the minimum 
energy configurations found are presented in Fig [5} Solutions were verified using a finite 
temperature annealing method. Particles are coloured according to their phases (see 
Fig [3]). 

We find that the point-particle approximation we have derived also favours configura¬ 
tions that form concentric ring-like structures. Furthermore, the qualitative forms of the 
energy minima predicted by the point-particle approximation are very close to the forms 
of the full numerical solutions found in Fig [lj although the results of the point-particle 
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Figure 3: Diagram illustrating the geodesics connecting pairs of particles on the disc, 
with the top pair being in phase (the maximally repulsive channel) and the bottom pair 
being out of phase (the maximally attractive channel). Particles are coloured according 
to their internal phase, using the bar on the right, with units given in multiples of n. 



Numerical - Morse 


Figure 4: Numerical and analytic approximations for the inter-soliton interaction. The 
upper curves are the results for baby Skyrmions in phase, and the lower curves are for 
baby Skyrmions out of phase. The blue curves are given by the Morse potential (21) with 
D = 0.76, p e = 0.73 and a = 1.13. Energies are given in units of 47T. 
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Figure 5: Minimal energy configurations for the point-particle approximation for topo¬ 
logical charges 1 < B < 20. Particles are coloured according to their internal phases (see 

Fig ||). 
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(a) B = 26 (b) B = 27 (c) £ = 28 



(d) B = 26 (e) 5 = 27 (f) B = 27* (g) B = 28 


Figure 6: The top row shows minimal energy configurations for the point-particle ap¬ 
proximation for topological charges 26 < B < 28. A popcorn transition to a three-layer 
solution is clearly seen at B = 27. This is confirmed by full numerical field calculations 
(below) which demonstrate a popcorn transition around B = 27, 28. At B = 27 two 
local minima with the same energy (to five significant figures) were found, indicated by 
the asterisk. 



(a) B = 53 (b) B = 54 



(c) B = 93 (d) B = 94 (e) B = 95 


Figure 7: Minimal energy configurations of the point-particle approximation demonstrat¬ 
ing popcorn transitions to four and five layers at charges B = 54 and B = 95 respectively. 
The lower charges pictured show how the inner-most ring of each configuration becomes 
deformed just before a pop. 
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charge B 

E/AttB 

form 

26 

1.8357 

{9, 17} 

27 

1.8546 

(9, 18} 

27* 

1.8546 

{1, 9, 17} 

28 

1.8723 

{1, 9, 18} 


Table 2: Energies per soliton and forms for AdS baby Skyrmions with topological charge 
26 < B < 28. Two local minima of the same energy (to five significant figures) were 
found at B = 27, as indicated by the asterisk. 


approximation are much more symmetric. This is to be expected since solitons in the full 
model have a finite size, and so can overlap and interact in more complicated ways. 

For charges B < 7 the approximation accurately predicts not only the ring structure 
of the configurations, but also the alignment of internal phases. For even and odd charges 
the particles have internal phase differences of 7r and n ± n/B respectively, as observed 
in the full model. At B = 8 the approximation gives the correct distribution of phases, 
although it does not reproduce the deformed ring structure. 

Furthermore, the point-particle approximation correctly captures the first popcorn 
transition at B = 9, and closely estimates the qualitative forms of the ring structures for 
all B < 20. The point-particle minima disagree with the full numerical results at charges 
B = 11, 13, 15, 16, 19 and 20. This may be due to the point-particle approximation 
not taking into consideration the radial forms of the B = 2 and 3 solitons, or may be a 
result of assuming the solitons can be approximated by particles with zero size. However, 
even when the exact forms do not agree, the difference is only by the position of a single 
particle. 

These results suggest that the point-particle approximation may be a useful tool 
in qualitatively estimating the forms taken by AdS baby Skyrmions for higher charges. 
Performing further numerical minimisations of (22) allows us to predict a second popcorn 
transition at charge B = 27. In fact, by using the predicted forms as a guide for choosing 
initial conditions, full numerical energy minimisations reveal the popcorn transition to 
three layers around charge B = 27, 28, although it is difficult to say exactly when the 
transition occurs due to the presence of two local minima with similar energies (shown 
in Fig |6j with energies given in Table [2]). 



(a) B = 200 (b) B = 250 

Figure 8: The minimal energy configuration of the point-particle approximation at 
charges B = 200 and 250. The ring structure near the origin is no longer visible, and 
may indicate the emergence of a lattice structure. 
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The point-particle approximation predicts a third and fourth popcorn transition at 
charges B = 54 and B = 95 (see Fig [7]). Since these charges are very large it would 
be difficult to perform an extensive search for the global minima in the full numerical 
model, even with the guidance provided by the point-particle approximation. However, 
the previous success of the model would indicate that the true popcorn transitions would, 
indeed, be near these points. 

Finally, investigation of the results predicted by the point-particle approximation for 
still higher charges may be able to provide clues as to the lattice structure preferred by the 
baby Skyrmions in the infinite charge limit. The energy minimum for charge B = 200 
can be seen in Fig [8j While clear rings can still be observed in the outermost layers, 
the centre of the configuration is heavily deformed and may indicate an emerging lattice 
structure. 


8 Conclusions 

We have investigated the static solitons and multi-solitons of the massless baby Skyrme 
model in a (2+l)-dimensional Anti de-Sitter spacetime. We have found that the spacetime 
curvature acts by adding an effective mass to the model which allows us to find static 
solutions to the equations of motion without a mass term. Solitons for topological charges 
1 < B < 3 were found to have a radially symmetric form, while higher-charge multi- 
solitons were found to form concentric ring-like solutions. 

As B increases, a series of transitions occur where the minimal energy solutions take 
the form of concentric rings with increasing numbers of layers, a phenomenon reminiscent 
of the baryonic popcorn transitions studied recently in the context of holographic dense 
QCD. In order to investigate these transitions further a point-particle approximation was 
derived which was able to qualitatively estimate the forms of minimal energy solutions 
for a wide range of topological charges, as well as accurately predict the charges at which 
further popcorn transitions occur. 

The point-particle approximation derived may indicate an emerging lattice structure 
for baby Skyrmions in AdS in the limit B —> oo. The minimum energy configuration 
found for B = 200 displayed clear rings towards the edge of the space, although near 
the origin the structure appeared significantly less ring-like, and seemed suggestive of 
an emergent lattice formation. Further investigation into this area would be required to 
make any further claims. 

The 0(3)-sigma model stabilised by a baby Skyrme term has previously been stud¬ 
ied in an AdS-like spacetime as a low-dimensional toy model of holographic QCD [Q], 
specifically the Sakai-Sugimoto model. It has been argued that a better low-dimensional 
analogue of this model would involve stabilising an 0(3)-sigma model using a vector me¬ 
son term, and such a model has been studied in a parameter regime where the two toy 
models are similar Pi. although investigation of a more interesting parameter regime 
proved difficult. It would therefore be interesting to study the vector meson model in 
pure AdS to see if a parameter regime could be found to give qualitatively different results 
to the baby Skyrme model. 

Finally, the natural extension to this paper is to study the full (3 + l)-dimensional 
Skyrme model in AdS. We have demonstrated that a multi-ring like structure exists in 
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two dimensions for AdS baby Skyrmions and this property may translate to the higher 
dimensional model. This could manifest as spherical multi-shells with polyhedral sym¬ 
metry groups and hence be approximated by multi-shell rational maps and it would be 
very interesting to study AdS Skyrmions to see if such configurations give lower energy 
solutions. 
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